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The fast-time wake transport and decay models require vertical profiles of crosswinds, 
potential temperature and the eddy dissipation rate as initial conditions. These inputs are 
normally obtained from various field sensors. In case of data-denied scenarios or 
operational use, these initial conditions can be provided by mesoscale model simulations. In 
this study, the vertical profiles of potential temperature from a mesoscale model were used 
as initial conditions for the fast-time wake models. The mesoscale model simulations were 
compared against available observations and the wake model predictions were compared 
with the Lidar measurements from three wake vortex field experiments. 
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Nomenclature 

vortex circulation (m 2 /s) 

initial vortex circulation (m 2 /s) 

initial vortex descent velocity (m/s) 

initial vortex pair separation (m) 

dimensional Brunt- Vaisala frequency (Y 1 ) 

dimensional eddy dissipation rate (m 2 Is 3 ) 

potential temperature ( K) 

temperature (°F) 

dew point temperature (°F) 

east- west velocity component (m/s) 

north-south velocity component (m/s) 


I. Introduction 

F AST-TIME wake transport and decay models are being developed and tested to provide solutions for safe and 
efficient aircraft wake spacing standards. They are also needed for safety and capacity benefit studies of 
advanced air traffic management concepts and procedures. Most of the fast-time models (e.g., Sarpkaya 2000; 
Sarpkaya et al. 2001; Holzapfel 2003; Proctor et al. 2006; Proctor 2009) require vertical profiles of temperature or 
potential temperature, crosswinds, and eddy dissipation rates (EDR) as initial conditions. These input weather 
conditions are usually obtained from various field sensors. For example, the cross winds can be obtained from the 
Lidar or an instrumented meteorological tower, the temperature from microwave radiometers, and the EDR from 
sonic anemometers. A vertical profile of EDR can be generated based on the similarity theory (Han et al. 2000) 
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provided that measurements are made at two different altitudes (e.g., from sonic anemometers mounted on a 
meteorological tower). Lidars can also be used to estimate the EDR (Pruis and Delisi 2011). Ideally, arrays of 
sensors could be deployed in the airport terminal area to provide spatially-varying information of different 
meteorological quantities. However, the cost of deploying a comprehensive array of sensors that can adequately 
cover the entire spatial extent of the terminal area is prohibitive and therefore measurements can typically be taken 
at just a few points in space. In addition, sensors cannot be used to forecast weather conditions several hours in 
advance. 

Weather forecasts for supporting the terminal area operations can be obtained with the help of numerical models 
and substantial progress has been made in mesoscale modeling over the last three decades (Ahmad and Proctor 
2012). The operational mesoscale weather models include the National Center for Atmospheric Research’s (NCAR) 
Weather Research and Forecasting (WRF) model (Klemp et al. 2007), Science Applications International 
Corporation’s Operational Multiscale Environment model with Grid Adaptivity (Bacon et al. 2000), the United 
States Navy’s Coupled Ocean/ Atmosphere Mesoscale Prediction System (Hodur 1997), and the National Centers for 
Environmental Prediction’s (NCEP) Rapid Update Cycle (RUC) developed by Benjamin et al. (2004). 

Several studies have been conducted in the past to evaluate mesoscale modeling for supporting terminal area 
operations (Kaplan et al. 1999; Charney et al. 2000; Kaplan et al. 2006; Freeh et al. 2007; Gerz et al. 2007; Ringley 
et al. 2007; Ahmad and Proctor 2012). These studies have primarily focused on quantifying errors in mesoscale 
model forecasts. In the current study, the WRF mesoscale model was used to generate initial conditions of potential 
temperature for the fast-time wake vortex transport and decay models. The simulations were conducted for three 
different wake vortex field experiments. The fast-time models were run using initial conditions from the WRF 
forecast and from observations. The fast-time model predictions were then compared with Fidar data. In the 
following sections, the results of fast-time wake vortex models and the WRF simulations are described in detail. 

II. Fast-Time Wake Transport and Decay Models 

Fast-time wake models are empirical or semi-empirical algorithms used for real-time predictions of wake 
transport and decay based on aircraft parameters and ambient weather conditions. The aircraft dependent 
parameters include the initial vortex descent velocity (V 0 ) and vortex pair separation distance ( b 0 ). The atmospheric 
initial conditions include vertical profiles of either temperature or potential temperature (0), EDR (s) and crosswind. 
The atmospheric parameters that affect wake decay are atmospheric stratification and turbulence. The first fast-time 
wake transport and decay model was developed by Greene (1986). Some of the fast-time models currently in use 
include: AVOSS (Aircraft Vortex Spacing System) Prediction Algorithm (APA), TASS (Terminal Area Simulation 
System) Derived Algorithms for Wake Prediction (TDAWP), Probabilistic 2-Phased (P2P) model and the 
Deterministic wake Vortex Model (DVM). The APA model (Robins and Delisi 2002) incorporates Sarpkaya’s 
model (Sarpkaya 2000; Sarpkaya et al. 2001) for out-of-ground-effect vortex decay and descent, and follows an 
approach similar to Corjon and Poinsot (1996) for in-ground-effect vortex transport. The TDAWP model (Proctor et 
al. 2006; Proctor and Hamilton 2009) has been developed from parametric studies using large eddy simulation of 
wake vortices. The APA and TDAWP models have been developed by the National Aeronautics and Space 
Administration (NASA), while the P2P model (Holzapfel 2003) and the DVM (de Visscher et al. 2010 ) have been 
developed by the Deutsches Zentrum ftir Fuft- und Raumfahrt (DFR) and Universite catholique de Fouvain (UCF), 
respectively. 


III. Mesoscale Modeling 

The WRF model is a mesoscale numerical weather prediction model that has been developed by NCAR with 
help from various partners in academia and federal research laboratories. It has governing equations for the time- 
dependent, three-dimensional, nonhydro static, fully compressible Navier-Stokes equations. The model grid 
structure is based on Arakawa C-grid staggering (Arakawa and Famb 1977). In the vertical, WRF uses a terrain- 
following pressure formulation, where the top of the domain is a constant pressure surface. The grid can be 
stretched in the vertical to provide higher resolution in the boundary layer. In the horizontal direction, one-way and 
two-way nesting options are available to provide high grid resolution in regions of interest. The horizontal and 
vertical advection terms can be discretized using schemes with second to sixth-order spatial accuracy. The time 
integration is done using an explicit third-order Runge-Kutta time marching scheme, with time-splitting of acoustic 
and gravity wave modes to maintain numerical stability. Choices for monotonic as well as positive -definite schemes 
are available for the advection of scalars. The latest version of WRF (Version 3.4) includes the option for a fifth- 
order Weighted Essentially Non-Oscillatory (WENO) scheme (Fiu et al. 1994) for the advection of microphysical 


2 

American Institute of Aeronautics and Astronautics 



scalars and turbulence kinetic energy (TKE). Extensive physics packages for modeling the surface layer 
interactions, boundary layer turbulence processes, and cumulus parameterizations have been implemented in the 
model. Different cloud microphysics schemes of varying complexity are available, as well as parameterizations for 
radiation physics. 


IV. Memphis Wake Vortex Field Experiment 


A. Memphis Field Experiment 

A comprehensive field experiment to measure wake vortices and the ambient meteorological conditions was 
conducted at the Memphis International Airport in Memphis, Tennessee from August 6 through August 29, 1995 
(Campbell, et al. 1997). The experiment was sponsored under NASA Langley Research Center’s Aircraft Vortex 
Spacing System (AVOSS) project (Hinton 1995; Perry et al. 1997). Several sensor systems were deployed to collect 
detailed meteorological data to study the effects of atmospheric conditions on the behavior of wake vortices. The 
sensors included radiosondes, sodars, a wind profiler, one 150 ft high meteorological tower, a Radio Acoustic 
Sounding System (RASS) and NASA Langley’s OV-10 research aircraft. The radiosondes were used to measure 
winds and temperature measurements (10s averages) at 50m vertical resolution. The OV-10 aircraft was flown at 
selected times and took measurements of temperature and winds at a sample rate of 10 Hz. Temperature (5 min 
averages) was measured using RASS every 30 min at 97 vertical levels from 145m to 1492m. The 150 |/i (45.7m) 
meteorological tower was equipped with a large array of sensor systems. Winds, temperature and moisture were 
measured from the tower at 5m, 10m, 20m, 30m, and 42m heights. Turbulence quantities (turbulence kinetic energy 
and eddy dissipation rate) were estimated from wind measurements at 5m and 40m heights. Rain rate, soil 
temperature, soil moisture, barometric pressure, and incoming and outgoing solar radiation also were measured by 
the sensors deployed on the meteorological tower. Standard meteorological data such as atmospheric pressure, 
temperature, moisture, cloud cover, and visibility, etc. were obtained from the National Weather Service’s Surface 
Aerodrome Observations (SAO) and the Automated Surface Observations System (ASOS). 

B. WRF Simulation Setup 

The WRF simulation domain consisted of an outermost domain bounded between 124.06°W and 55.93°W in 
longitude and 18.08°N and 49.14°N in latitude, with a horizontal mesh resolution of 36 km. Two higher resolution 
nests were defined within the outermost domain. Domain 2 was bounded between 111.54°W and 68.45°W in 
longitude and 24.84°N and 43.86°N in latitude with a horizontal mesh resolution of 12km. The innermost domain 
was bounded between 99.74°W and 79.97°W in longitude and 30.66°N and 39.29°N in latitude with a horizontal 
mesh resolution of 4 km. High-resolution (1km) terrain and land use datasets were used for the innermost domain 
while lower resolution terrain was used in the coarser nests. A stretched grid with 60 levels was used in the vertical 
to provide high-resolution grid spacing in the boundary layer. The WRF computational domain for Memphis is 
shown in Figure 1 . 

The fifth-order upwind-biased scheme was used in the horizontal and the third-order upwind-biased scheme was 
used in the vertical within a three-stage Runge-Kutta explicit time-marching scheme. The monotonic scheme option 
was used for the transport of both the microphysical scalars and the TKE. The Rapid Radiative Transfer Model 
(RRTM) longwave radiation scheme (Mlawer et al. 1997) and the Goddard shortwave scheme (Chou and Suarez 
1999) were used to parameterize the effects of both longwave and shortwave radiative transfer in the atmosphere. 
The Thompson scheme (Thompson et al. 2004) was used for the microphysics, and the cumulus parameterization 
schemes were switched off for the simulations. Mellor-Yamada-Janjic turbulence scheme was used to parameterize 
subgrid scale turbulent processes. The simulation was initialized for August 15, 1995 at 1800UTC using the North 
American Regional Reanalysis (Mesinger et al. 2006) data from the National Climate Data Center (NCDC) and was 
run for a 30 hr forecast (spin-up of 6 hr). Data from the high-resolution innermost nest were output every 15 min. 

C. Comparison of WRF Predicted Temperature and Winds with Observations 

The WRF simulation results were evaluated against available observations. Figure 2 shows the comparison of 
simulated surface temperatures and winds with the ASOS data and the SAVPAK sensors mounted on the met tower. 
Forecast errors were quantified in terms of root-mean- square error ( rmse ), mean absolute error ( mae ), and mean 
bias. These errors for the twenty-four hour period at Memphis are listed in Table 1. Overall, the WRF model 
performed well and was able to simulate the main features of the diurnal cycle for August 16, 1995. The surface 
temperature, dew point temperature and wind predictions from the model were in good agreement with observations 
at Memphis, both qualitatively as well as quantitatively. 
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Figure 1. The WRF computational domain with higher resolution nests. The innermost nest had a resolution 
of 4km and was centered at Memphis, Tennessee. 



time (utc) 


<D 

u 

-p 

nJ 

U 

<1) 

CU 

W 


d) 

V 

-P 

c 

■rl 

o 

o* 

* 


d) 

T> 



u (sav) 


v(sav) 


u (wrf ) 


v(wrf) 

■ 

u ( obs ) 

■ 

v ( obs ) 


to 



time (utc) 


Figure 2. Comparison of WRF simulations with observations. WRF simulation results are denoted by solid 
lines. Observed temperatures at 2m for Memphis are shown as squares. Observed winds at 10 m for 
Memphis are depicted by squares. The SAVPAK 5m data is given by small diamonds in the figure. The 
SAVPAK data frequency was every one minute. 


Table 1: WRF Surface Forecast Errors for Memphis, Tennessee (August 16, 1995) 



T(°F) 

T d (°F) 

u(m/s) 

v(m/s) 

case 

rmse 

mae 

bias 

rmse 

mae 

bias 

rmse 

mae 

bias 

rmse 

mae 

bias 

1995-08-16 

1.40 

1.15 

0.81 

1.37 

1.02 

-0.66 

1.68 

1.32 

0.74 

0.88 

0.74 

-0.43 


4 

American Institute of Aeronautics and Astronautics 


v (m/s) 



D. Fast-Time Wake Transport and Decay Model Results 

There were a total of 46 landings for August 16, 1995, in the AVOSS database. Wake transport and decay 
predictions were obtained from three fast-time models (APA3.2, APA3.4, and TDP2.1) for each of the events using 
known aircraft parameters and inputs of crosswind, temperature, and turbulence. Two sets of simulations were 
conducted for each fast-time model: 1) the vertical profile potential temperature from observations was used, and 2) 
the vertical profile of the WRF-predicted potential temperature was used. The potential temperature profiles in the 
AVOSS database were generated by fusing data from several sensors. A comparison of the average Brunt-Vaisala 
frequency from the observations and the WRF forecast is shown in Figure 3. 

The fast-time model predictions were then evaluated using Continuous Wave (CW) Lidar observations. The 
results were quantified in terms of root-mean- square error, and are given in Table 2. The circulation error was 
normalized by the initial circulation strength and the errors in lateral and vertical transport were normalized by the 
initial vortex separation, b 0 . All three models had relatively low errors for circulation and lateral transport. 

A comparison of the TDP2.1 wake predictions with CW Lidar data is shown in Figures 4-5. The generating 
aircraft in this case was a DC 10 and the time corresponds to the nighttime stable boundary layer. The TDP2.1 
simulations (initialized by observations and the WRF data) of circulation decay and vortex descent compare well 
with the Lidar measurements. 
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Figure 3. Comparison of the average Brunt-Vaisala frequency, N (s' 1 ) for each vertical profile obtained from 
the WRF simulation and the observations. Memphis - August 16, 1995. 


Table 2: Fast-Time Models Evaluation (Mean RMSE) - Memphis: August 16, 1995 (46 Cases) 



Circulation 
(normalized by T 0 ) 

Lateral Transport 
(normalized by b 0 ) 

Altitude 

(normalized by b 0 ) 


ICs from OBS 

0 from WRF 

ICs from OBS 

0 from WRF 

ICs from OBS 

0 from WRF 

TDP2.1 

0.211 

0.210 

0.706 

0.704 

0.488 

0.484 

APA3.2 

0.165 

0.167 

0.645 

0.644 

0.506 

0.501 

APA3.4 

0.162 

0.162 

0.641 

0.640 

0.504 

0.503 
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Figure 4. Comparison of TDP2.1 predicted wake vortex time histories with the observed CW Lidar data. 
The generating aircraft was a DC10 in this case. The left panel shows TDP2.1 simulation initialized 

by observed potential temperature. The right panel shows TDP2.1 simulation initialized by WRF predicted 
potential temperature. Memphis - August 16, 1995. 




Figure 5. Comparison of TDP2.1 predicted wake vortex time histories with the observed CW Lidar data. 
The generating aircraft was a DC10 in this case. The left panel shows TDP2.1 simulation initialized 

by observed potential temperature. The right panel shows TDP2.1 simulation initialized by WRF predicted 
potential temperature. Memphis - August 16, 1995. 
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V. Dallas-Fort Worth Wake Vortex Field Experiment 

A. Dallas-Fort Worth Field Experiment 

In 1997, the AVOSS (Hinton 1995; Perry et al. 1997) setup at Memphis International Airport was relocated to 
Dallas-Fort Worth International Airport. The primary objective of this field deployment was to conduct a 
comprehensive test of various AVOSS procedures established at Memphis (Dasey et al. 1998; Joseph et al. 1999). 
Wake measurements were obtained from NASA’s 2 pm pulsed Lidar, the Lincoln Labs’ 10.6pm Continuous Wave 
Lidar and an array of anemometers (windline). One of the objectives was to collect meteorological data for the 
validation of mesoscale models and therefore an impressive array of meteorological sensors was deployed to 
quantify the atmospheric state in detail. Two meteorological towers were installed during this field experiment. 
Each meteorological tower was equipped with various sensor systems for the measurement of temperature, relative 
humidity, wind speed, and wind direction at different heights. Turbulence measurements were taken with sonic 
anemometers which measured the three components of wind at 10 Hz. Several sensors were mounted at the base of 
each tower to measure the atmospheric pressure, rainfall, and solar radiation, etc. In addition to the two 
meteorological towers, a wind profiler, a RASS and a sodar were also deployed. Special upper air soundings were 
taken six times a day from five different sites surrounding the Dallas-Fort Worth International Airport. Standard 
meteorological data such as atmospheric pressure, temperature, moisture, cloud cover, and visibility, etc. were 
obtained from the National Weather Service’s SAO and the ASOS data. 

B. WRF Simulation Setup 

Three levels of grid nesting were used in the simulation in order to achieve acceptable grid resolution while 
maintaining computational efficiency. The WRF simulation domain consisted of an outermost domain bounded 
between 129.16°W and 64.91°W in longitude and 13.44°N and 50.10°N in latitude with a horizontal mesh 
resolution of 36 km. The innermost high resolution domain was centered at Dallas-Fort Worth, Texas and bounded 
between 102.73°W and 90.67°W in longitude and 29.26°N and 36.92°N in latitude with a horizontal mesh 
resolution of 4 km. High-resolution (1km) terrain and land use datasets were used for the innermost domain while 
lower resolution terrain was used in the coarser outer nest. A stretched grid with 60 levels was used in the vertical to 
provide high-resolution grid spacing in the boundary layer. The WRF computational domain is shown in Figure 6. 

The Mellor-Yamada-Janjic turbulence scheme (Mellor and Yamada 1978; Janjic 1994) was used to parameterize 
sub-grid scale processes and higher-order schemes were used for computing the convective fluxes (fifth-order 
upwind-biased in the horizontal and third-order upwind-biased in the vertical). A fifth-order WENO scheme with 
positive definite limiter was used to transport microphysical scalars as well as turbulence kinetic energy. The option 
of WENO schemes is available in the latest distribution of the WRF model (Version 3.4). The simulation was 
initialized using the North American Regional Reanalysis (Mesinger et al. 2006) data from the NCDC at 0000UTC 
on October 3, 1997 and was run for a 24 hr forecast. Data from the high-resolution innermost nest were output every 
15 min. 

C. Comparison of WRF Predicted Temperature and Winds with Observations 

Four previous studies have focused on the meteorological aspects of the AVOSS deployment at Dallas-Fort 
Worth (Kaplan et al. 1999; Charney et al. 2000; Kaplan et al. 2006; Ringley et al. 2007). Kaplan et al. (1999) 
discussed the formation of the Great Plains low level nocturnal jet and the ability of mesoscale models to accurately 
simulate the phenomenon. Charney et al. (2000) and Ringley et al. (2007) evaluated the accuracy of turbulence 
kinetic energy and EDR prediction from mesoscale models. In Kaplan et al. (2006) surface temperature and dew 
point temperature errors (~1°C) were estimated from surface maps for one data point only on two different days 
(September 18, 1997 1200UTC and September 19, 1997 1200UTC). An error analysis for the complete diurnal 
cycle using surface observations was not reported in any of these studies. In this study, the WRF simulation was run 
specifically to provide initial conditions of potential temperature for the fast-time models and therefore, the WRF 
simulation results were first evaluated against available observations. Figure 7 shows the comparison of simulated 
surface temperatures and winds with the ASOS data. Forecast errors were quantified in terms of root-mean- square 
error ( rmse ), mean absolute error ( mae ), and mean bias. These errors for the twenty-four hour period at Dallas-Fort 
Worth are listed in Table 3. Overall, the WRF model preformed reasonably well for October 3, 1997. The errors for 
the Dallas-Fort Worth simulation were larger for temperature and dew point temperature (compared to the Memphis 
simulation) but similar for the u and v components of the wind. 
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Figure 6. The WRF computational domain with higher resolution nests. The innermost nest had a resolution 
of 4 km and was centered at Dallas-Fort Worth, Texas. 




Figure 7. Comparison of WRF simulations with observations. WRF simulation results are denoted by solid 
lines. Observed temperatures at 2m and observed winds at 10 m for Dallas-Fort Worth are shown as squares. 


Table 3: WRF Surface Forecast Errors for Dallas-Fort Worth, Texas (October 3, 1997) 



T(°F) 

T d (°F) 

u(m/s) 

v(m/s) 

case 

rmse 

bias 

mae 

rmse 

bias 

mae 

rmse 

bias 

mae 

rmse 

bias 

mae 

1997-10-03 

4.11 

3.87 

3.87 

4.72 

-3.30 

3.70 

1.98 

1.43 

1.58 

0.93 

-0.39 

0.76 
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D. Fast-Time Wake Transport and Decay Model Results 

There were a total of 42 landings for October 3, 1997 for Dallas-Fort Worth in the AVOSS database. Wake 
transport and decay predictions were obtained from three fast-time models (APA3.2, APA3.4, and TDP2.1) for each 
of the events using known aircraft parameters and inputs of crosswind, temperature, and turbulence. The potential 
temperature profiles in the AVOSS database were generated by fusing data from several sensors. Two sets of 
simulations were conducted for each fast-time model: 1) the vertical profile potential temperature from observations 
was used, and 2) the vertical profile of the WRF -predicted potential temperature was used. A comparison of the 
average Brunt- Vaisala frequency from the observations and the WRF forecast is shown in Figure 8. 

The fast-time model predictions were then evaluated using Continuous Wave (CW) Lidar observations. The 
results were quantified in terms of root-mean- square error, and are given in Table 4. The events on October 3, 1997 
were mostly in near-ground effect (NGE). In the proximity of ground, stratification plays a secondary role in 
circulation decay of wakes. This is reflected in almost identical results obtained using the observed potential 
temperature and WRF predicted potential temperature (Table 4). 

A comparison of the TDP2.1 wake predictions with CW Lidar data is shown in Figures 9-10. The generating 
aircraft in this case was a B752. The TDP2.1 predictions (initialized by observations and the WRF data) of 
circulation decay and vortex altitude compare well with the Lidar measurements. 


0.03 


0.025 



Figure 8. Comparison of the average Brunt-Vaisala frequency, N (s' 1 ) for each vertical profile obtained from 
the WRF simulation and the observations. Dallas-Fort Worth - October 3, 1997. 


Table 4: Fast-Time Models Evaluation (Mean RMSE) - Dallas-Fort Worth: October 3, 1997 (42 Cases) 



Circulation 
(normalized by T 0 ) 

Lateral Transport 
(normalized by b 0 ) 

Altitude 

(normalized by b 0 ) 


ICs from OBS 

0 from WRF 

ICs from OBS 

0 from WRF 

ICs from OBS 

0 from WRF 

TDP2.1 

0.233 

0.233 

0.792 

0.792 

0.278 

0.278 

APA3.2 

0.245 

0.245 

0.728 

0.728 

0.267 

0.266 

APA3.4 

0.229 

0.229 

0.728 

0.728 

0.270 

0.270 
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Figure 9. Comparison of TDP2.1 predicted wake vortex time histories with the observed CW Lidar data. 
The generating aircraft was a B752 ^ ^ in this case. The left panel shows TDP2.1 simulation initialized 
by observed potential temperature. The right panel shows TDP2.1 simulation initialized by WRF predicted 
potential temperature. Dallas-Fort Worth - October 3, 1997. 




Figure 10. Comparison of TDP2.1 predicted wake vortex time histories with the observed CW Lidar data. 
The generating aircraft was a B752 - 4 in this case. The left panel shows TDP2.1 simulation initialized 

by observed potential temperature. The right panel shows TDP2.1 simulation initialized by WRF predicted 
potential temperature. Dallas-Fort Worth - October 3, 1997. 
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VI. Denver Wake Vortex Field Experiment 


A. Denver Field Experiment 

The Denver field experiment was conducted by NASA during the month of September 2003. Although, the 
primary objective for this deployment was to evaluate wake measurements using acoustic sensors (Dougherty et al. 
2004), wake trajectories and circulation were also measured using pulsed and continuous wave Lidars. In addition, 
temperature profiles were measured with the MTP5 sensor (microwave radiometer), and wind profiles from the 
pulsed Lidar measurements. During the Denver experiment, one minute average winds were measured by propeller 
anemometers that were mounted on a meteorological tower. The altitudes of the sensors were at 7m, 14.6m, and 
32.3m AGL. In addition, wind measurements at a rate of 10 Hz were obtained with an ultrasonic anemometer that 
was mounted on a 7.3m pole. Compared to the AVOSS deployments at Memphis (1995) and Dallas-Fort Worth 
(1997), the meteorological data collection at Denver was not as comprehensive. The generation of a vertical profile 
of eddy dissipation rate requires two data points at two different heights (Han 2000). In Denver 2003 experiment, 
there was only one sonic anemometer and therefore the profiles of eddy dissipation rates were obtained from the 
Lidar data using spatial structure functions (Pruis and Delisi 2011). 

B. WRF Simulation Setup 

The WRF simulation domain consisted of an outermost domain bounded between 130.25°W and 79.65 °W in 
longitude and 24.48°N and 53.24°N in latitude with a horizontal mesh resolution of 12km. The inner high resolution 
domain was centered at Denver, Colorado and bounded between 115.71°W and 93.52°W in longitude and 32.77°N 
and 46.48°N in latitude with a horizontal mesh resolution of 4 km. High-resolution (1km) terrain and land use 
datasets were used for the innermost domain while lower resolution terrain was used in the coarser outer nest. A 
stretched grid with 90 levels was used in the vertical to provide high-resolution grid spacing in the boundary layer. 
The WRF computational domain is shown in Figure 1 1 . 

The Mellor-Yamada-Janjic turbulence scheme (Mellor and Yamada 1978; Janjic 1994) was used to parameterize 
sub-grid scale processes and higher-order schemes were used for computing the convective fluxes (fifth-order 
upwind-biased in the horizontal and third-order upwind-biased in the vertical). Fifth-order WENO scheme with 
positive definite limiter was used to transport microphysical scalars as well as turbulence kinetic energy. The 
simulation was initialized using the North American Regional Reanalysis (Mesinger et al. 2006) data from the 
NCDC at 0000UTC on September 19, 2003 and was run for a 24 hr forecast. Data from the high-resolution inner 
nest were output every 15 min. 

C. Comparison of WRF Predicted Temperature and Winds with Observations 

Several simulations were conducted for Denver (September 14 - September 20) in order to minimize the forecast 
errors. First, a set of sensitivity runs was conducted which included spinup time, microphysics scheme, radiation 
scheme, number of vertical levels, and horizontal mesh resolution. Spinup times of 3 hr, 6 hr, and 12 hr degraded the 
solution and therefore the simulation was initialized at 0000UTC. The choice of microphysics scheme other than 
Thompson did not improve the forecast and in some cases degraded the solution. The choice of radiation schemes 
and horizontal mesh resolution did not improve the model forecast. 

The WRF Data Assimilation Package (WRFDA) was also evaluated for Denver simulations. Data from the 
National Oceanic and Atmospheric Administration’s (NOAA) Meteorological Assimilation Data Ingest System 
(MADIS) and surface and upper air observations from the NCAR archives were used in this study. The MADIS 
dataset included metars, satellite observations, mesonet data, raobs, Aircraft Communications Addressing and 
Reporting System (ACARS) data, and NOAA profilers. WRF options of Objective Analysis, Four-Dimensional 
Data Assimilation (Observational Nudging) and 3D Variational Analysis (Barker 2004) were tested. For the 
September 19, 2003 simulation, the four-dimensional data assimilation option resulted in an improvement in the 
forecast. 

Figure 12 shows the comparison of simulated surface temperatures and winds with the ASOS data. The WRF 
data assimilation was used in this simulation. The forecast errors for the twenty-four hour period at Denver are 
listed in Table 5. The baseline simulation listed in Table 5 was initialized using gridded data only and the 
baseline+on simulation used four-dimensional data assimilation with observational nudging. The errors in 
temperature and especially dew point temperature were reduced by using data assimilation. However, the error in v 
component of wind increased slightly. Figure 12 shows that much of the error in temperature forecast was due to 
WRF’s failure to capture the rapid decrease in temperature during the nighttime. 
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Figure 11. The WRF computational domain with higher resolution nest. The innermost nest had a 
resolution of 4 km and was centered at Denver, Colorado. 
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Figure 12. Comparison of WRF simulations with observations. WRF simulation results are denoted by 
solid lines. Observed temperatures at 2m and observed winds at 10 m for Denver are shown as squares. 


Table 5: WRF Surface Forecast Errors for Denver, Colorado (September 19, 2003) 



T(°F) 

T d (°F) 

u(m/s) 

v(m/s) 

case 

rmse 

bias 

mae 

rmse 

bias 

mae 

rmse 

bias 

mae 

rmse 

bias 

mae 

baseline 

4.18 

2.47 

3.10 

4.46 

-2.73 

3.98 

1.91 

-0.35 

1.20 

1.40 

-0.55 

1.14 

baseline + on 

3.87 

2.15 

3.21 

3.69 

3.05 

3.19 

2.03 

-0.01 

1.50 

2.08 

-0.64 

1.54 
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D. Fast-Time Wake Transport and Decay Model Results 

A total of 82 landings for September 19, 2003 at Denver were used in this evaluation. The vertical profiles of 
temperature at Denver were measured by the MTP5 sensor, which is a microwave radiometer. 

Two sets of simulations were conducted for each fast-time model: 1) the vertical profile potential temperature 
from observations was used, and 2) the vertical profile of the WRF -predicted potential temperature was used. A 
comparison of the average Brunt- Vaisala frequency from the observations and the WRF forecast is shown in Figure 
13. The fast-time model predictions were then evaluated using CTI Pulsed Lidar observations. The results were 
quantified in terms of root-mean- square error, and are given in Table 6. All three models had relatively low errors 
for circulation and lateral transport. 

A comparison of the TDP2.1 wake predictions with CTI Pulsed Lidar data is shown in Figures 14-15. The 
generating aircraft in this case was an MD83. The TDP2.1 predictions (initialized by observations and the WRF 
data) of circulation decay and vortex descent compare well with the Lidar measurements. 
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Figure 13. Comparison of the average Brunt- Vaisala frequency, N (s' 1 ) for each vertical profile obtained 
from the WRF simulation and the observations. Denver - September 19, 2003. 


Table 6: Fast-Time Models Evaluation (Mean RMSE) - Denver: September 19, 2003 (82 Cases) 


Model 

Circulation 
(normalized by T 0 ) 

Lateral Transport 
(normalized by b 0 ) 

Altitude 

(normalized by b 0 ) 


ICs from OBS 

0 from WRF 

ICs from OBS 

0 from WRF 

ICs from OBS 

0 from WRF 

TDP2.1 

0.211 

0.228 

0.616 

0.626 

0.874 

0.868 

APA3.2 

0.218 

0.209 

0.620 

0.621 

0.933 

0.913 

APA3.4 

0.184 

0.195 

0.613 

0.615 

0.926 

0.930 
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Figure 14. Comparison of TDP2.1 predicted wake vortex time histories with the observed CTI Pulsed 

Lidar data. The generating aircraft was an MD83 in this case. The left panel shows TDP2.1 

simulation initialized by observed potential temperature. The right panel shows TDP2.1 simulation 
initialized by WRF predicted potential temperature. Denver - September 19, 2003. 




Figure 15. Comparison of TDP2.1 predicted wake vortex time histories with the observed CTI Pulsed 

Lidar data. The generating aircraft was an MD83 — in this case. The left panel shows TDP2.1 
simulation initialized by observed potential temperature. The right panel shows TDP2.1 simulation 
initialized by WRF predicted potential temperature. Denver - September 19, 2003. 
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VII. Summary 

In this study, the state-of-the-art in mesoscale modeling was evaluated for supporting terminal area operations in 
general and wake vortex applications in particular. The National Center for Atmospheric Research’s WRF 
mesoscale model was used to conduct simulations for three different geographical areas (Memphis, Tennessee; 
Dallas -Fort Worth, Texas; and Denver, Colorado). The WRF simulation data were evaluated against available 
observations with favorable results. The errors in surface temperature, dew point temperature, and wind speed and 
direction were relatively small for all three locations. The WRF forecasted stratification also compared well with 
observations - with the exception of Denver. In the case of Denver, WRF was unable to simulate the rapid cooling 
at surface during the nighttime and the level of stratification in the nighttime stable boundary layer. 

Two sets of fast-time model simulations were conducted: 1) initialization using observed potential temperature 
profiles, and, 2) potential temperature profiles from WRF forecasts. The fast-time model predictions were compared 
against the Lidar measurements. In general, all three fast-time models (APA3.2, APA3.4, and TDP2.1) gave similar 
results with relatively low errors. The errors in fast-time models predictions were also low when initialized by WRF 
predicted potential temperature and comparable to errors in predictions initialized by observations. This study 
demonstrates the feasibility of using high-resolution mesoscale simulations for fast-time wake vortex modeling 
applications. 


Appendix A 

The model inputs specified in the simulations were based on the WRF sensitivity study conducted by Ahmad and 
Proctor (2012). The main input parameters used in the simulations are listed in this appendix. 


Table A.l: Model Options used in WRF Simulations 


Model Option 

Memphis (MEM) 

Dallas/Fort Worth (DFW) 

Denver (DEN) 

WRF Model Version 

3.2.1 

3.4.1 

3.4.1 

Initialization 

1995-08-15 (1800UTC) 

1997-10-03 (0000UTC) 

2003-09-19 (0000UTC) 

Forecast Hours 

30 hr 

24 hr 

24 hr 

Ax = Ay (Innermost nest) 

4 km 

4km 

4km 

Vertical Levels 

60 

60 

90 

Microphysics Scheme 

Thompson 

Thompson 

Thompson 

PBL Scheme 

Mellor-Yamada-Janjic 

Mellor-Yamada-Janjic 

Mellor-Yamada-Janjic 

Scalar Advection 

Monotonic 

5 ,h -order WENO 

5 th -order WENO 

Initialization Data 

NARR 

NARR 

NARR 

Data Assimilation 

No 

No 

Yes 
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